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We present a showcase for deriving bounds on the neutrino masses from laboratory experiments 
and cosmological observations. We compare the frequentist and Bayesian bounds on the effective 
electron neutrino mass mp which the KATRIN neutrino mass experiment is expected to obtain, using 
both an analytical likelihood function and Monte Carlo simulations of KATRIN. Assuming a uniform 
prior in mp, we find that a null result yields an upper bound of about 0.17 eV at 90% confidence 
in the Bayesian analysis, to be compared with the frequentist KATRIN reference value of 0.20 eV. 
This is a significant difference when judged relative to the systematic and statistical uncertainties 
of the experiment. On the other hand, an input mp = 0.35 eV, which is the KATRIN 5a detection 
threshold, would be detected at virtually the same level. Finally, we combine the simulated KATRIN 
results with cosmological data in the form of present (post-WMAP) and future (simulated Planck) 
observations. If an input of mp = 0.2 eV is assumed in our simulations, KATRIN alone excludes a 
zero neutrino mass at 2.2a. Adding Planck data increases the probability of detection to a median 
2.7a. The analysis highlights the importance of combining cosmological and laboratory data on an 
equal footing. 

PACS numbers: 14.60.Pq, 02.50.Tt, 95.85.Ry 



> 

m 

On 
O 

o 



X 



One of the most important recent discoveries in physics 
is that neutrinos have non-zero masses. Since neutrinos 
are explicitly massless in the Standard Model, massive 
neutrinos necessarily involve new physics and may point 
the way for other developments. Therefore there is much 
need for improved knowledge of neutrino masses from 
both laboratory experiments and cosmological observa- 
tions. 

The fact that neutrinos are massive is demonstrated by 
atmospheric 0, 0, solar W i, IS IS 0, reactor @, i, [l0| 
and accelerator fill. [lH . Il3| neutrino oscillation exper- 
iments. Additionally, the masses are constrained from 
above by tritium beta decay [3, UH, neutrinoless dou- 
ble beta decay 0, H3, EI El , and fr om bounds on the 
amount of hot dark matter in the Universe (see below 
for references). Since the oscillation experiments are only 
sensitive to square mass differences, the absolute neutrino 
mass scale has not been determined, apart from the dis- 
puted observation of neutrinoless double beta decay in 
the Heidelberg-Moscow experiment [IS [U [22j. How- 
ever, recent cosmological observations of the Cosmic Mi- 
crowave Background [H, l2i. l25l, Igol \27\ , the Large Scale 
Structure of the Universe , the Baryon Acoustic 

Oscillations in the matter power spectrum [3011 . Ty pe la 



Supernovae [3l|, l32|, |33[ , and Lyman-a clouds [34j show 
substantial evidence for the 'standard' ACDM cosmolog- 
ical model. In that model case, the upper bound on the 
sum of neutrino masses ^2 m v is driven well into the sub- 



cV region [27|, |35j, |3fl l37|, |3J, |3J, |40| . If so, it seems un- 
likely that a detection of the electron neutrino mass mp 
will be made at KATRIN [4l[, the next generation tri- 
tium beta decay experiment which has a discovery poten- 
tial of 0.35 eV. However, there are two important issues 
with the cosmological data. Firstly, the results depend on 
which cosmological model is used, and different models 
contain different numbers of parameters. Adding more 
parameters will in general relax the bounds on any sin- 
gle parameter. Secondly, it should be noted that most 
of the predicting power relies on Lyman-a data and it is 
extremely hard to model these systems |4(J. This makes 
it challenging to assess the systematic effects associated 
with this measurement. Without Lyman-a, current lim- 
its are in the region 0.6-0.7 eV for ^ m v . Hence a cleaner 
probe such as laboratory experiments, galaxy surveys or 
weak lensing surveys is very desirable. It should be noted 
that the next generation of large-scale surveys should be 
able to access the range down to 0.2 eV for the sum of 
the neutrino masses [42|. 

Comparison of various neutrino mass bounds is com- 
plicated by the fact that there are several different, but 
interdependent, observables involved (see, e.g., Fogli et 
al. [II), but also by the use of different statistical meth- 
ods which potentially leads to inconsistencies. The cos- 
mology community usually adheres to Bayesian inference 
while experimental neutrino physicists tend to employ 
frequentist methods in oscillation experiments and abso- 
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lute mass searches. 

In this paper we address these complications by com- 
paring frcqucntist and Baycsian neutrino mass bounds in 
relation to the near-future KATRIN experiment and by 
combining the estimated KATRIN results with present 
and future cosmological observations in a consistent man- 
ner. First, we introduce the KATRIN experiment and 
then we compare the frequentist and Baycsian confidence 
intervals obtained from an approximate KATRIN likeli- 
hood function. Then we simulate KATRIN using Monte 
Carlo methods and compare a Bayesian analysis of our 
simulations with the frequentist analysis by the KATRIN 
collaboration. Finally, we predict the future neutrino 
mass sensitivity by combining our results with present 
(post-WMAP) and future (Planck) cosmological obser- 
vations in a consistent Bayesian framework. 



I. THE KATRIN EXPERIMENT 

The Karlsruhe Tritium Neutrino experiment KATRIN 
is a next generation single beta decay experiment investi- 
gating the kinematics of the tritium beta spectrum near 
its kinematical endpoint [4l[ . It will measure the effective 
electron neutrino mass through the nuclear beta decay of 
tritium 



H 



5 He 4 



(1) 



with an endpoint energy Eq ~ 18.6 keV. Any non-zero 
(effective) mass of the electron antincutrino would lead 
to an end of the beta spectrum lower than the endpoint 
energy Eq as well as to a distortion of the spectrum in 
the region near the endpoint. 

KATRIN investigates the /3-electron energy with an 
electrostatic high-pass filter which measures the inte- 
grated spectrum of /3-electrons near the endpoint of the 
spectrum. The lower bound of the filter is determined by 
a retarding potential U and the number of electrons in a 
given potential energy bin is [43j 



N(qU) = A 



dE- 



dN 
dE 



(E,E o ,m 2 )f ies (qU,E) + N b , 



(2) 

where f rcs is the experimental response function, Eq is 
the endpoint energy of the /3-spectrum, N b is the back- 
ground, and A is the signal amplitude of the spectrum. 
The differential electron spectrum close to the endpoint 
is, approximately, 



dN 
dE 



C{Eq-E)^{Eq-EY 



(3) 



The observable is therefore the square of the effective 
mass. The aim of the KATRIN experiment is a sensitivity 
of m/3 < 0.2 eV at 90% confidence in case of a null result 
or a 5tr discovery potential for nip > 0.35 eV (43J. It has 
been shown 44[ that non-SM contributions to the beta 
decay such as right-handed currents have little impact on 
the neutrino mass measurement. 



KATRIN is actually being built at the site of 
Forschungszentrum Karlsruhe, having major components 
such as the main spectrometer vessel with its diameter of 
10m and a length of 24m already installed. Commission- 
ing of the 70m-long full setup comprising the windowless 
gaseous tritium source, the electron transport section, a 
tandem spectrometer and a semiconductor detector array 
is expected for 2010. To reach the expected sensitivity, a 
measuring time of 3 full-beam-years is anticipated, start- 
ing end of 2010. 

Here, our interest in KATRIN is mainly its expected 
statistical performance. The benchmark is the frequen- 
tist Monte Carlo analysis by the KATRIN collaboration 
as given in the Design Report [43|: The outcome is an 
expected distribution of occurrence which is very close to 



Gaussian in mi with a width of cr sia: t.( m g) 



0.018 eV^ 



This is combined with the estimate of the total system- 
atic uncertainty, cr sys .( m2 ) < 0.017 eV 2 , which leads to 

T 2 _i_ ^2 
slat. 



the total uncertainty of atot.( m ^) 



a" Y' 2 



0.025 eV 2 . Taking the Gaussian shape of the frequency 
distribution, the 90% confidence sensitivity on mp is then 
derived via 



L(90%) - Vl.64(7 to t., 



(4) 



and the detection threshold is nip = [5o"tot. { m2 p)] 1 ^ 2 = 
0.35 eV. 



II. BAYESIAN VS. FREQUENTIST 
CONFIDENCE INTERVALS 

Next we discuss the frequentist and Bayesian ap- 
proaches to parameter inference and calculation of upper 
and lower bounds. 

Both frequentist confidence intervals and Bayesian 
credible intervals are constructed from a likelihood func- 
tion L(x, s) which is a function of the data x and the 
parameter to be inferred s (For simplicity, we consider 
only one parameter in this discussion). A frequentist 
confidence interval is a statement about frequency of oc- 
currence. The interval [si,S2] is a member of the set of 
all intervals that fulfil the condition that the probability 



P{s € [si,s 2 ]) = 1 -a, 



(5) 



for all s, in particular the unknown true value. Here, I— a 
is the stated confidence level of the interval. The values 
Si and s 2 arc determined by the likelihood function. This 
defining property known as 'coverage' ensures that, in the 
limit of a large ensemble of intervals, a fraction 1 — a of 
the ensemble will contain the true value. 

In the unified approach of Feldman and Cousins 
[45l ] . intervals are constructed from a likelihood order- 
ing method which involves two steps. Firstly, for a fixed 
choice of s a physically allowed value of x is admitted to 
an 'acceptance region' if the likelihood at that value is 
greater than at any other value not included in the inter- 
val already. In the continuous case, this means that the 
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acceptance region defined by [0:1,0:2] will satisfy 



L(x, s)dx = 1 



(6) 



so that, additionally, the likelihood function L is every- 
where greater inside the region than outside. Secondly, 
this procedure is repeated for every s on a finely spaced 
grid. Having measured the best-fit value x , the con- 
fidence interval is the union of all the values of s that 
include xq in their acceptance region. 

On the other hand, the Bayesian credible interval is 
a statement about degree of belief. It attaches a proba- 
bility to a range of values, indicating quantitatively how 
much it is believed that the actual value is within that 
range. Following Bayes' theorem, the posterior probabil- 
ity density in the true value of the parameter s is 



P{s\x) = 



L(x\s)ir(s) 



(7) 



where L{x\s) is the likelihood function, which is now 
viewed as the conditional probability of the data given 
s, and 7r(s) is the assumed prior probability density of 
the parameter. The normalizing factor P{x), also known 
as the Bayesian evidence, is simply the integral of the 
numerator. The credible interval [si,S2] is calculated by 
integrating the posterior probability density so that 



1 



P(x) 



L(x\s)ir(s)ds = 1 



(8) 



The framework can straightforwardly be generalized 
to multiple parameters, including nuisance parameters 
which are easy to handle. One can just integrate 
(marginalize) the multidimensional posterior distribution 
over the nuisance parameters to obtain the posterior in 
the parameter (s) of interest. The prior is a means of in- 
corporating known information about s, which is then 
updated by the measurements in the form of the likeli- 
hood. One can think of at least three bases for the choice 
of prior: a known experimental result, a theoretical anal- 
ysis, or 'ignorance' about the value of the parameter. 



III. COMPARISON OF INTERVALS FOR 
KATRIN 

Now we study the frequentist and Bayesian intervals 
constructed from the predicted KATRIN likelihood func- 
tion for a vanishing mass. As discussed, the observable 
is the square of the electron neutrino mass mjg and for 
the purposes of fitting it is necessary to also consider 



negative values of this parameter, i.e. 



< 0. The 



realness is then imposed afterwards. Therefore we as- 
sume a likelihood which is Gaussian in m \ and has width 
.(ml) = 0.025 eV 2 , the estimated total error of KA- 
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TRIN ^ 

the range ±5er = ±0.125 eV 2 around zero. We calcu- 
late the frequentist confidence intervals in the unified 



|43| . We consider varying best-fit values f in 



— Feldman-Cousins 
■ - Flat prior in m 2 

— Flat prior in m 
Flat prior in log{m) 




0.00 
Best-fit m 2 , eV ! 



FIG. 1: Upper and lower bounds at 90% confidence as a func- 
tion of the best-fit value calculated in the frequentist Feldman- 
Cousins approach, and in the Bayesian approach with three 
different priors. The interpretation of a null result, m 2 , f = 
is very sensitive to the choice of analysis, while the confidence 
interval for a frequentist 5<r detection at mj j = 0.125 eV 2 is 
virtually the same in all cases. A Gaussian likelihood function 
in m 2 a is assumed. Note that, for the purpose of fitting, the 
observable m 2 } can take on both positive and negative values. 



approach of Feldman & Cousins [45j but we state the 
bounds in m^, the parameter of interest, rather than in 



We also calculate Bayesian credible intervals for 
three different priors which are constant in m 2 }, m^, and 
log 77^3, respectively. The first two priors are bounded 
from below at zero mass and the logarithmic prior is 
bounded at 7 meV, suggestive of lower limits from os- 
cillation data. All the priors are bounded from above at 
1 eV. 

There is a certain freedom in calculating the Bayesian 
intervals, since an interval containing a specific proba- 
bility mass is not uniquely determined without further 
constraints. We choose to use the shortest possible in- 
terval, integrated in the metric in which the chosen prior 
is constant. However, we stress that this choice is not 
unique. 

The resulting upper and lower bounds at 90% con- 
fidence are shown in Fig. [T] The intervals, regardless 
of method and assumed prior, approximately converge 
when the best-fit value is physical and far from the 
boundary, i.e. for to 2 f > 0.125 eV 2 = ha. The be- 



haviour for to 2 f < varies considerably; in particular 
the Feldman-Cousins upper bound drops faster than any 
Bayesian upper bound when going from positive to neg- 
ative m 2 , f . Also, the Bayesian upper bounds, except for 
the logarithmic prior, are more conservative for very neg- 
ative best-fits. The behaviour of the logarithmic prior is 
distinct due to the larger prior probability density at- 
tached to small values of 771,3. This causes the logarith- 
mic prior to yield the tightest upper bounds for small 
to 2 f . At first glance the logarithmic prior, or Jeffreys' 
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prior [46[ as it is also known, may appear artificial but 
it has the nice property that it treats all powers of the 
parameter equally. Also, as shown by Jaynes [43], it is 
the prior which contains the least information about a 
scale parameter, i.e. an 'ignorance' prior. 

Finally, we note that the question of whether a nonzero 
neutrino mass is detected can also be discussed from a 
model selection perspective. In the frequentist frame- 
work this is the rejection of the null hypothesis that the 
mass is zero, while in the Bayesian sense one may com- 
pare the Bayesian evidence for models with and without 
the neutrino mass 14a. 



IV. BAYESIAN BOUNDS FROM SIMULATIONS 
OF KATRIN 

In order to understand better the statistical pcrfo- 
mance of KATRIN from a Bayesian point of view we 
perform an independent Monte Carlo simulation of the 
experiment which is analyzed in the Bayesian framework 
and compared with reference results. 

The model is essentially given by eqs. ([2]) and j3]). 



There are four input parameters: 



l 8' 



the endpoint en- 



ergy of the beta spectrum Eq, the background per bin 
A r b, which is assumed to be independent of U, and the 
signal amplitude of the spectrum, A. We adjust A to re- 
produce the statistical error of the frequentist KATRIN 
simulations. We also recover the approximately Gaussian 
shape of the likelihood function. 

We assume priors which are uniform in each of the 
nuisance parameters and in mp. For each realization of 
the data we calculate the likelihood function on a grid 
in the four dimensional parameter space which extends 
more than five standard deviations to each side of the in- 
put values. The grid contains 50 points along each of the 
three nuisance parameters and 130 points along the nip 
direction in the range [—0.1, 1] cV 2 and it is constructed 
so that it is regular in mp. From this tabulated likelihood 
function we marginalize over the three nuisance parame- 
ters and convolve the result with a Gaussian with width 
a S ya(m 2 p) = 0.017 eV 2 to model the maximum expected 
systematic error |43| . The result is the posterior prob- 
ability density of mp from which the Bayesian credible 
interval can be calculated as in section HT1 

We have checked that the resolution of the grid is 
sufficient by calculating the total probability contained 
within a given interval for ten datasets resolved on our 
grid and on a reference grid with 100 points for each nui- 
sance parameter. The results differ by 0(1O~ 5 ). 

We consider two cases with input masses mp = and 
mp = 0.35 eV, respectively. For the zero input we cal- 
culate the upper bound and for the nonzero input the 
central value of the confidence interval, both at 90% con- 
fidence. The resulting distributions are shown in Figs. [2] 
and O respectively. Also shown are the expected upper 
bounds/central values (dotted lines) obtained in section 
HT1 as well as the frequentist expectation. 



Upper bound on neutrino mass 
Input mass m.=0 eV 



Bayesian 
Feldman-Cousins 




0.10 0.15 
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FIG. 2: Distribution of 90% upper bounds on mp for an input 
mp = 0, obtained in the Bayesian analysis assuming a uniform 
prior in mp. As expected, the distribution is scattered around 
the upper limit from Fig. [1] for a Gaussian likelihood with a 
uniform prior in mp, which is significantly lower than the 
frequentist expectation. Out of 1000 simulated datasets, 67 
did not include zero at 90% confidence. 
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Bayesian 
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0.30 0.35 0.40 

Center of 90% credible interval, eV 



FIG. 3: Distribution of central values of 90% confidence inter- 
vals for an input mp — 0.35 eV. As in Fig. [5] the distribution 
is scattered around the central value for the corresponding in- 
terval for a uniform prior in Fig. [1] but in this case there is no 
discrepancy with the frequentist expectation for the central 
value. 



The mean upper bound of the distribution in Fig. [2] is 
0.17 eV so it is clear that the Bayesian analysis is more 
optimistic than the frequentist which yields the bound 
0.20 eV. This is in agreement with the analysis in section 
HIl in particular with the upper bound of 0.175 eV ob- 
tained for a uniform prior in rap and best-fit m 2 f = 0, as 
shown in Fig. [T] This difference between the frequentist 
and Bayesian upper bound is large since, if the frequen- 
tist analysis were to yield an upper bound of 0.17 eV, 
the total error of <Jtot.{ m2 p) would need to be reduced by 
40%. 
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FIG. 4: The level at which mp = is just excluded, given 
in terms of the equivalent number of cr's for the input mp = 
0.35 eV. The left vertical line is the conventional 5<r detection 
threshold and the 63% of the datasets to the right of this line 
exclude mp = at more than 5a. The right vertical line is the 
expected exclusion level in the Bayesian analysis, assuming a 
uniform prior in mp and a Gaussian likelihood in m 2 p. 51% 
of the datasets are to the right of this line. 



For the input mp = 0.35 eV in Fig.[3j the central value 
distribution is centered on the input and the Bayesian 
and frequcntist analyses yield the same result. Again, 
this is in agreement with the intervals in Fig. [T] since the 
input mass is much greater than zero. The crucial ques- 
tion for the discovery potential, however, is at which con- 
fidence level the simulated datasets just exclude mp = 0, 
in terms of the equivalent number of cr's. This is shown 
in Fig. 3] in terms of the equivalent number of cr's. Also 
indicated is the conventional 5tr discovery limit. It can 
be seen that a majority of the datasets exclude mp = 
at a greater confidence level than 5er. The reason for 
this is mainly the choice of prior; a uniform prior in m 2 ^ 
instead would only show random scatter around the 5tr 
line. When assuming a likelihood function which is Gaus- 
sian in mp and adopting a uniform prior in mp, rap = 
is excluded at a probability corresponding to 5.2a. This 
value is very close to the median of the distribution in 
Fig.H 

From a frequcntist point of view Fig. 2] indicates that 
the Bayesian credible intervals with uniform prior in mp 
are more likely to undercover. This means that less than 
the fraction 1 — a of an ensemble of intervals will contain 
the true value of the parameter — at least at this extreme 
confidence level. 

In summary, the Bayesian analysis carried out here is 
more optimistic than the frequcntist in this respect in the 
sense that both the sensitivity and discovery potential are 
lower. 



FIG. 5: Relation between mp and ^ m v for normal and in- 
verted hierarchies, and the naive relation ^m v = 3mp. For 
mp > 0.1 eV there is virtually no differences between the 
relations. 

V. COMBINATION WITH COSMOLOGY 

Having analyzed the differences between mass bounds 
derived by the different statistical methods, we proceed 
by combining our KATRIN simulation with present and 
future cosmological observations in the Bayesian frame- 
work. 

The observable in tritium beta decay mp and the cos- 
mological observable ^ m„, are related through the neu- 
trino mixing angles and square mass differences which 
are probed by oscillation experiments. The mixing an- 
gles form a parametrization of the PMNS matrix which 
relates the weak interaction neutrino eigenstates \v a ) to 
the mass eigenstates \vi) through \v a ) = U*i \ vi) [49j . 
We simplify the analysis by fixing the oscillation data 
[1, H, El at best-fit values. In particular, we assume 
sinf9 13 = which determines mp = (mi + W2)/2, and 
the masses can then be written in terms of observables 
as 

m\ = m 2 p - ^Am 2 12 , (9a) 

ml = m 2 p + ^Am 2 12 , (9b) 

ml = m 2 p + \^Am\ 2 ± Am 2 3 . (9c) 

The only freedom left is the hierarchy, i.e. the sign of 
Am 2 3 in (|9c| . The best-fit square mass differences arc 

Am\ 2 = 7.9 x 1CT 5 eV 2 , (10a) 
Am| 3 = 2.6 x 10" 3 eV 2 . (10b) 

The hierarchy determines if the mass m-$ is greater or 
lesser than the pair m\ and m2- Relations between mp 
and m v are shown in Fig. O 

In principle, the uncertainties on the neutrino mixing 
angles and square mass differences should be taken into 
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Sum of neutrino masses Lm, eV Equivalent number of o's 



FIG. 6: The one-dimensional marginalized posteriors for 
52 m v as obtained in in the simulation of Planck, computed 
with an input mass of input Yl m v = 0.6 eV, and the 
WMAP analyses. These posteriors yield upper bounds of 
Y^rriu < 1.1 eV (Planck low), Jm„ < 0.81 eV (Planck 
high), < 1.8 eV (WMAP only), and ]Tm„ < 0.66 

eV (WMAP+LSS+SNIa), all at 95% confidence. 



FIG. 7: Exclusion levels of mp — given in the equivalent 
number of ct's for an input mp = 0.2 eV to the KATRIN simu- 
lation, with and without cosmological data. The distribution 
for KATRIN alone differs from Fig.|4]only through the input. 
The extra data from cosmology clearly improve the detection 
threshold, although for this input mass a 4a detection is still 
unlikely. 



account, but for the purpose of this showcase we deem 
the best-fit values sufficient. It was found in ref. [13] that 
the uncertainty in the m v relation is very small, 

except for the smallest allowable masses. 

We combine simulated KATRIN posteriors with 
present cosmological data in the form of the WMAP 3- 
year data to which data from large-scale structure sur- 
veys and type la supernovae are added. Additionally, 
we also combine the KATRIN results with a simulation 
of the Planck satellite, the next cosmic microwave back- 
ground satellite. 



A. Bounds from the CMB 

We include present limits on m„ from cosmology 
in the form of the analysis of the WMAP 3-year data, as 
performed by the WMAP team [27[. We use their MCMC 
chains, available at [53], for a ACDM cosmology de- 
scribed by the six minimal parameters plus the neutrino 
energy density. We use the chains for the largest com- 
bined dataset, which includes WMAP and other CMB 
data as well as large scale structure and supernovae mea- 
surements. It was shown in rcf. [13] that an identified 
problem with the WMAP likelihood code does not af- 
fect neutrino mass bounds from combined datasets. On 
its own, this data set yields an upper bound on the 
sum of the neutrino masses of J2 m f < 0-66 eV at 95% 
confidence. This strong bound is mainly caused by the 
large-scale structure and small-scale CMB data included; 
the WMAP data alone yield the much weaker constraint 
Y J m v < 1.8 eV, also at 95% confidence [13, H3] ■ 

Additionally, we have obtained forecasts from the 



Planck satellite by using the technique described in Ab- 
dalla & Rawlings [42| . This consists mainly in assuming 
a fiducial model and calculating the expected rejection 
probability of a model around this fiducial value. We 
have assumed conservative values for the sensitivity by 
taking only one usable science channel with 8 detectors, 
a noise effective temperature of 120 (iK^/s, an angular 
resolution of 10 arcseconds, and 60% sky coverage in a 
one year survey. We have also considered a more opti- 
mistic forecast which would cover around 70% of the sky 
in a two year survey. In both cases the input fiducial mass 
was taken as m„ = 0.6 eV and the other cosmologi- 
cal parameters considered were: O m = 0.3, Sit, = 0.04, 
n s = 0.95, ft A = 0.7, n run = 0.00, cr 8 = 0.8, r = 0.09, 
and we assumed a flat Universe. The resulting posteriors 
yield upper bounds of m v < 1.1 eV (conservative) and 

m„ < 0.81 eV (optimistic) at 95% confidence. 

Both the WMAP+LSS+SNIa posterior and the simu- 
lated Planck posteriors, marginalized over all parameters 
except ^Jm„, are shown in figure [5] The Planck simu- 
lation yields higher upper bounds on mv since it is 
based on the CMB alone, and because the input to the 
simulation is already close to the upper bound derived 
from WMAP+LSS+SN. 



B. Combined KATRIN and CMB 

Now we combine our simulation of KATRIN with the 
cosmological observations discussed above. 

We simulate a new set of 1000 KATRIN experiments 
with an input mass of mp = 0.2 eV, the expected KA- 
TRIN 90% upper bound on a null result. Using the rela- 
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tions in ((9]), we convert each simulated m^-posterior into 
a posterior in ^ m v , assuming normal hierarchy for def- 
initcncss. For simplicity, we do not propagate the small 
errors from the oscillation data. The simulated KATRIN 
posteriors are combined with each of the cosmology pos- 
teriors individually and for each of these we calculate the 
confidence level at which a 'vanishing' neutrino mass is 
excluded. There is a technical point here since converting 
the KATRIN observable using (|9|) implies that neutrinos 
are massive. In other words, whereas it is possible to 
examine a KATRIN posterior and a cosmology poste- 
rior in isolation and determine whether there is evidence 
for a nonzero neutrino mass in each posterior, the two 
can only be combined under the assumption that both 
m/3 and ^2 "ii/ are nonzero. Therefore we search for the 
confidence level at which the lower bound of the inter- 
val just reaches the lower limit of ^ m v and regard this 
as the detection threshold (for that posterior) of a mas- 
sive neutrino. Otherwise the detection of the nonzero 
mass is from the oscillation data rather than the KA- 
TRIN+cosmology posterior. 

The procedure yields distributions of the exclusion 
level of a 'vanishing' mass for KATRIN alone and for 
the two types of combined posteriors, which are shown 
in Fig. [71 For the KATRIN posteriors alone, the median 
detection level is 2.2a and it is clear that a large frac- 
tion of the simulated experiments do not exclude mp = 
even at 2a. On the other hand, most experiments do 
exclude mp = at 90% confidence, in agreement with 
expectations based on the analysis in section ITT1 

The addition of the cosmological data improves the 
dctcctability of a nonzero mass. The combination with 
WMAP+LSS+SNIa yields a median detection level of 
2.6a while the combination the Planck yields 2Aa (con- 
servative) and 2.7a (optimistic). More importantly, the 
widths of the distributions become smaller so that many 
more combined posteriors exclude a vanishing mass at 
more than 2a. On the other hand, there are no detec- 
tions at more than 4<r which is a consequence of the fact 
that neither of the three datasets are able to provide 
significant detections on their own. The combination 
of the high specs-Planck simulation with KATRIN ac- 
tually performs slightly better than WMAP+LSS+SNIa 
and KATRIN, which is likely due to the fact that the 
WMAP+LSS+SNIa posterior is peaked at a much lower 
value of m v than the input 0.6 eV used in the simula- 
tions. In summary, the combination of KATRIN and cos- 
mological data may well extend the discovery potential of 
the experiments into regions not accessible by either type 



of experiment alone. We note that additional large-scale 
structure data could improve the detection potential fur- 
ther. 



VI. CONCLUSIONS 

We have presented a showcase for estimating the neu- 
trino masses by combining data from laboratory and cos- 
mological experiments. Specifically, we combine simu- 
lated results from KATRIN with cosmological probes. 
A key issue in this context is the consistent application 
of statistical methods to all data. We analyze the dif- 
fering results obtained by the frequentist and Bayesian 
approaches, taking simulated KATRIN likelihoods as an 
example. Then wc combine the KATRIN results with 
both the present cosmological constraints after WMAP 
and simulated future data from the Planck satellite. 

The comparison of frequentist and Bayesian methods 
shows important differences for KATRIN, where the like- 
lihood is likely peaked close to the physical boundary 
mp = 0. While the KATRIN reference upper bound to a 
null result is 0.2 eV at 90% confidence, the Bayesian re- 
sult is 0.17 eV, assuming a uniform prior, and even lower 
if a logarithmic prior is assumed. As explained in section 
IIV1 this is a sizable difference judged against the system- 
atic and statistical uncertainties. Far away from the zero 
mass boundary, the different approaches converge. 

The combined constraints from laboratory and cosmo- 
logical data are investigated by assuming an electron neu- 
trino mass of 0.2 eV. In this case KATRIN alone will 
provide a 2.2ct detection, based on 10 3 simulations of 
KATRIN. The 10 3 exclusion levels span a rather wide 
range, from zero detection to about Aa. The addition 
of CMB data increases the detection level somewhat to 
2.6cr (WMAP+LSS+SNIa) or 2Aa - 2.7er (Planck) but, 
importantly, narrows the range of exclusion levels so that 
in most cases a > 2a detection is found. 

Our analysis highlights the impact of the choice of sta- 
tistical analysis, and it demonstrates the application of 
Bayesian methods to combined datasets. 
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